The few-body problem for trapped bosons with large scattering length 
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We calculate energy levels of two and three bosons trapped in a harmonic oscillator potential with 
oscillator length flosc- The atoms are assumed to interact through a short-range potential with a 
scattering length a, and the short-distance behavior of the three-body wave function is characterized 
by a parameter 9. For large positive a/aosc, the energies of states which, in the absence of the trap, 
correspond to three free atoms approach values independent of a and 6. For other states the 9 
dependence of the energy is strong, but the energy is independent of a for |a/aoac| ^ 1- 
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By exploiting the properties of Feshbach resonances, 
it has become possible to tune atomic scattering lengths 
to essentially any value, either repulsive or attractive 0. 
This has led to renewed interest in the properties of sys- 
tems with large scattering length a, a problem which 
has also been studied extensively in the context of nu- 
clear physics, because the scattering lengths for nucleon- 
nucleon interactions are large in magnitude compared 
with the spatial scale tq of the interparticle interaction. 

A novel possibility for exploring experimentally the 
properties of few-body systems is to trap atoms in a pe- 
riodic potential created by standing-wave light beams, a 
so-called optical lattice. For example, in the experiments 
described in Ref. §, up to 200,000 *'^Rb atoms were 
distributed over more than 150, 000 lattice sites, thereby 
creating a large number of few-body systems with up 
to 2.5 atoms each on average. This provides one mo- 
tivation for theoretical studies of few-body problems in 
traps. A second motivation is to understand theoretically 
the properties of bulk matter when the magnitude of the 
scattering length is large compared with the interparticle 
spacing n~^/^, where n is the particle density. In this 
regime, the usual low density approximations, which are 
valid only if n|ap <C 1, fail. Recently, one has sought a 
universal description, with n\a\'^ as the only parameter, 
of the many-body problem for Fermi |^ and Bose Q sys- 
tems interacting via short-range potentials (ro <C n^^^^), 
but with n|a|^ ~ 1. Studies of three-body recombination 
in ultracold Bose gases have, however, revealed that the 
rate depends on an additional three-body parameter , 
and in a recent work, it has been suggested that the equa- 
tion of state of a Bose gas in the limit n^^^^ ^ \a\ may 
depend on the three-body parameter Study of the 
few-body system in a trap is therefore of interest in or- 
der to obtain insights into the properties of bulk matter. 

In this Letter we calculate the energy of systems with 
two or three identical bosons of mass m in a harmonic 
trapping potential Vtrap = mu;^r^/2, where uj is the fre- 
quency of oscillations of a particle about the minimum 
of the potential. The length associated with the zero- 
point motion of a particle in the trap is Gobc — \/h/muj. 
The bosons are assumed to interact via a short-range 



potential. To take this into account, we impose the 
boundary condition dhi{rij'^)/drij\r-^^o — — 1/a- Here 
\I'(ri, . . . Tat) is the wave function of the iV-body system, 
Ti being the coordinates of the particles, and = r — r^. 

We first consider the two-body problem. Solving the 
Schrodinger equation for two bosons in a harmonic os- 
cillator potential, one finds that the eigenenergies E, ex- 
cluding the center of mass contribution, are given by ||^ 
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For a — the s-wave states have energies E = 'hu){2v + 
3/2), (z^ — 0, 1, ... ) at which the gamma function in the 
denominator diverges, while for a/flosc ±oo, the en- 
ergies are E = huj{2i/ + 1/2). Hence, when the bosons 
are confined to a volume smaller than |ap, the scattering 
length becomes unimportant. Interestingly, the energy 
does not increase without limit as a becomes large, as 
one would predict from the standard low-density result 
for the energy. 

We now investigate the analogous problem for three 
trapped bosons, and adopt the method of Ref. The 
wave function of three identical bosons may be decom- 
posed in the usual Faddeev fashion as 

vl/(ri, r2, ra) = i^(x, y) + V-Cx', y') + ^^(x", y"). (2) 

Here the Jacobi coordinates x = ri2, and y = — 
(ri + r2)/2 have been used. If interactions between 
atoms are purely s-wave, ^ is a function of the hy- 
perradius p = (x^/2 -I- 2y^/3)^/^ and the hyperangle 
a = arctan(V3a;/2y) only. Variables related to the orig- 
inal ones by cyclic permutations of the particles are de- 
noted by primes and double primes. 

The wave function ijj may be represented exactly by an 
infinite sum 
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where (j)i{p,a) are eigensolutions to the hyperangular 
part of the problem, which for a short-range potential 
reduce to 



1 



20 




-5 
-10 
-15 



' 1 




_ Efimov limit 


XThree free atoms 


1 


1 IVIolecule + atom 
1 1 , 1 



-2 2 4 

Logio(p/|a|) 

FIG. 1. Eigenvalue \i{p/\a\) of the hyperangular equation 
for a > (solid) and a < (dashed). 



(j)i{p,a) = \jN{p) sin y/\i{a-TT/2) 



(4) 



The eigenvalues \i are determined by the boundary con- 
dition as a — > 0, i.e. when two of the bosons approach 
each other, and are given by the equation 
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The hyperradial wave functions fi are the solutions of 
an infinite number of coupled differential equations 
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The effective potential for the hyperangular motion 
therefore depends on the eigenvalue associated with 
the hyperangular degrees of freedom. The coupling ma- 
trices are given by the kinetic-energy operator acting on 
the hyperangular wave functions: 
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where <^i{p,a,a' ,a") = 0j(/3, a) -|- 0i(p, a') -I- ^^(p, a")- 
An approximation which has worked well in a number of 
applications |^ is to neglect the coupling terms P and 
Q. This is an adiabatic approximation in the sense that 
the hyperangular functions <i>i(/ci, a, a', a") are assumed 
to vary much more slowly with hyperradius than with 
the hyperangles. The method can be refined by including 
coupling between a finite subset of states. 

In Fig. |l| we plot \i{p/\a\) calculated from Eq. (||) for 
some adiabatic states. For both a > and a < we show 



the state approaching Ai = 4 for large p/|a|. This would, 
in the absence of the trap, correspond at large p to three 
free atoms, and it is the state of greatest interest in the 
present work. Also shown for a > is the result for the 
state which would correspond to a diatomic molecule and 
an atom. 

Let us now examine the form of the adiabatic hyper- 
spherical wave functions. The functions of interest are 
ones which remain finite as p ^ 0. The effective potential 
for small p/\a\ varies as (Ai(0) — 1/4)^"^. Consequently 
for Ai(0) > the hyperradial wave function fi tends to 
zero as pV^^^'W+i/s j^j. p/|a|, and it does not de- 

pend on the nature of three-body correlations when no 
two atoms are separated by more than the typical range 
of the interaction. The situation is quite different for 
Ai(0) < 0, since the effective potential is then attractive 
at small p. The two linearly independent solutions oscil- 
late rapidly as p — > 0, but do not diverge, and the general 
solution for small p may be written 



fip) \/psin[|so| ln(p/aosc) + 0], 



(8) 



where sq = Ao(0) ~ 1. 00624? and is a three-body 
parameter analogous to the phase shift for the two-body 
problem. 

We turn now to the energy eigenvalues. In Fig. ^ we 
exhibit numerical results for the solution of Eq. ^ , where 
the only states included are those shown in Fig. |l]. For 
a > 0, coupling between the two hyperangular functions 
has been included. Since this is weak, it is still possible 
to separate states that for large p correspond to three free 
atoms from states that correspond to a molecule and an 
atom. We show the energy only of the states that for 
large p correspond to three free atoms, since these are 
the states of interest here. 
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FIG. 2. Three-body energies as a function of a/ao 
0(mod7r) = (crosses) and ^(modTr) = 7r/2 (diamonds). For 
comparison the solid line shows the energy obtained when 
three-body correlations are neglected. 
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For small |a/aosc| the energy of the lowest harmonic 
oscillator state can be expanded as 



6 



(9) 



A lower bound to the energy may be found by omit- 
ting Qii{p) in Eq. (^). An upper bound may be ob- 
tained by evaluating the energy for the trial wave func- 
tion fi^^^i, where is the solution of Eq. (||) in the 
absence of the coupling terms Including the same 

states as before, we find —0.55 < c < —0.25. For com- 
parison, if three-body correlations are neglected, the cor- 
relation energy (the difference between the total energy 
and that of the corresponding state in the absence of 
interactions) is three times the correlation energy for 
two bodies, as obtained from Eq. (||). The resulting ap- 
proximation for the energy has the same term linear in 
d/cLosc, but the quadratic term has a different coefficient, 
c = 6(1 — ln2)/7r ~ 0.59. Hence, three-body effects are 
to leading order proportional to (a/aosc)^- 

With increasing |a|/aosc the dependence of the energy 
on 9 grows. For a < this is due to the sensitivity 
of the adiabatic state to the inner boundary condition. 
The adiabatic state for a > is insensitive to 9, and 
the 9 dependence arises due to coupling to states hav- 
ing different hyperangular functions. Energies for = 
and 9 = tt/2 are compared in Fig. ||, and the difference 
grows with a/aasc until a/flosc ~ 1, where it is ~ 5% 
of the correlation energy. At larger values of a/uosc the 
9 dependence decreases. This is an effect of the confin- 
ing potential. The hyperradial wave function /(p) varies 
on a scale Oosc- As shown above, the hyperangular 
wave function <i>i does, however, depend on the hyperra- 
dius through the combination p/a. Hence the couplings 
P and Q Eq. (0) are suppressed by factors Oosc/a and 
(aosc/a)^ respectively when a > Qosc- 

The 9 dependence of the energy is shown in more detail 
in Fig. P for a > 0. Physical conditions only determine 
0(mod7r), but to exhibit the structure we have continued 
the energy curves over a larger range of 9. Hence, given 
values of a/aosc and 9, E{9), E{9+'k), E{9+2'k) etc. along 
the same curve represent the energies of different phys- 
ical states. States corresponding mainly to free atoms 
at large p are represented by the nearly ^-independent 
segments of the curves. In the adiabatic approximation, 
i.e. neglecting all couplings, these states have energies 
shown by the dotted lines. Energies of other states, i.e. 
states associated with molecular bonding, are given by 
the dashed lines. The difference between the solid and 
dashed/dotted hues is the effect of the couphngs between 
the two hyperangular functions. The couplings are most 
important when a ~ Oosc, and are almost entirely ab- 
sent for < a < O.lflosc and a > lOflosc- The results 
for a = lOaosc indicate that for a ^ Oosc the couplings 
vanish and our solution becomes exact. 




FIG. 3. Energy of 9 for (from below) a/oosc = 0.1, 0.5, 1, 
2, and 10. Dashed and dotted lines show energies in the adia- 
batic approximation, i.e. neglecting couplings between hyper- 
angular functions. has been extended over several branches 
TT in order to simultaneously show several physical states (see 
text). 

In the limit a/aosc ^ ico Eq- @) may be solved ex- 
actly. The energy tends to a constant value, which for 
the lowest state with a given Ai(0) > is E = Tiuj{l + 
^ Ai(0)). For the state which at large p corresponds to 
free atoms, and for a > 0, one has Ai(0) ~ 19.94 and 
E ~ 'bAlhio. The corresponding energy if only two-body 
correlations are included \s E — Gtiuj, and therefore the 
correlation energy is reduced by approximately 18% by 
three-body effects. For Ai(0) < the energies are func- 
tions of the boundary condition for p — > 0, i.e. of 9, and 
are given by the solutions of the equation 
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For E > this equation gives harmonic oscillator 
states shifted by the interaction, with spacings En+i — 
En — 2fiu), while for i? < it gives bound three-body 
states, so-called Efimov states with energies scaling like 
En+i/En = exp(27r/|so|) — 515 |ll| perturbed by the 
oscillator potential. For example 9 = Q gives ener- 
gies E ~ -291649?k^, -566^^, -0.85?iw, 2.27nuj, 
4.48?iw, .... 

Let us now examine the implications of our results for 
the properties of bulk matter. On physical grounds, we 
would expect the energy per particle in trapped few-body 
systems to be similar to that for bulk matter if the oscil- 
lator length is replaced by the typical particle separation. 
For |a| 3> ftosc, the energy per particle of two- and three- 
body states for which Ai(0) > is ~ ?iw greater than in 
the absence of interactions. Since the density of particles 
is ri ~ l/a^gj., this corresponds to an energy per particle 
~ (T? jrn)v?l'^ . There is no dependence on 9 for these 
states. The state corresponding to unbound atoms in a 
trap falls into this class for positive a. This suggests that 
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the energy per particle of bulk matter has the universal 
dependence ^ Cn^/^, where C is a constant that does not 
depend on the three-body parameter 9. In particular, it 
does not have density-dependent contributions that are 
periodic in Inn, contrary to the possibility suggested in 
Ref. There are other states, such as those associated 
with molecular states bound by the interatomic potential, 
which do exhibit a strong dependence on 6. These lat- 
ter states are the final states in recombination processes, 
where diatomic molecules are formed from free atoms, 
and consequently recombination rates depend strongly 
on 9. 

The expansion Eq. of the energy of three trapped 
bosons, may be compared to the low-density expansion 
for the energy of N bosons in a homogenous gas |jl^ 
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For comparison, our result Eq. (^) for the energy of three 
trapped bosons with |a/aosc| ^ 1, rnay be rewritten as 
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Here the average density, weighted by the particle 
number, is defined hy n = J n{r)dN{r) / J dN{r) — 
J n^{r)dr/N, where dN{r) = n{r)dr, and therefore for 
a homogeneous condensate h = n. The three-body and 
many-body results agree for the term linear in a, which 
usually is the only term retained in, e.g. the Gross- 
Pitaevskii equation for condensates. Beyond the linear 
term the two expansions differ. The next term in Eq. 
( |rT| ) is of order (na^)^/^ relative to the first one, while 
in Eq. (|l^ ) it is proportional to a/oosc- If J; Oosc/i 
the most important correction to the linear dependence 
of energy on a will be the second term in Eq. ([l^), which 
arises from the trapping potential. 

To determine when the three-body parameter 9 be- 
comes important, we note that, in the absence of the 
trap, the S-matrix element between the adiabatic state 
corresponding to free atoms and that corresponding to 
a diatomic molecule and an atom is proportional to 
(fca)^, where k is the wave number of the colliding atoms. 
Therefore the recombination rate is proportional to a*, 
in agreement with earlier work Q . The constant of pro- 
portionality depends on 9, because of the 9 dependence 
of the final state. The same result is valid in a trap with 
lose 3> |a|, where due to the zero point motion k is of or- 
der l/flosc- Therefore, the first term in the small a/aosc 
expansion of the energy [Eq. (^] that depends explicitly 
on the coupling between different adiabatic states (and 
therefore on 9) is of order (a/aosc)"*- In the expansion Eq. 
( p^ ) this corresponds to the cubic term, since there one 
power of a/flosc has been factored out. To this order the 



result agrees with the results for the energy density of the 
homogeneous Bose gas in Ref. m. As we saw above, for 
a ^ aosc the couplings P and Q between hyperangular 
functions are reduced, which will tend to suppress recom- 
bination rates. 

In conclusion, we have found that there are states of 
the trapped three-boson system with very weak depen- 
dence of the energies on the three-body parameter 9. 
This dependence vanishes completely in the limit of an in- 
finitely large scattering length. The fact that the energy 
of these states is independent of the three-body param- 
eter in the limit of large scattering lengths gives some 
grounds to hope that the energy of trapped systems with 
a larger number of bodies is insensitive to the parameters 
describing the four- and higher-body problems, and that 
similar conclusions apply to the energy of bulk matter. 
Topics for future study include higher-body correlations, 
and the lifetime of a condensate with a large scattering 
length. 

We are grateful to Eric Braaten and Oleg I. Kartavtsev 
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